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STRONG SOLUTIONS TO THE COMPRESSIBLE 
NAVIER-STOKES-VLASOV-FOKKER-PLANCK EQUATIONS: 
GLOBAL EXISTENCE NEAR THE EQUILIBRIUM AND LARGE 

TIME BEHAVIOR 


FUCAI LI, YANMIN MU’, AND DEHUA WANG 


Abstract. A kinetic-fiuid model describing the evolutions of disperse two- 
phase flows is considered. The model consists of the Vlasov-Fokker-Planck 
equation for the particles (disperse phase) coupled with the compressible Navier- 
Stokes equations for the fluid (fluid phase) through the friction force. The fric¬ 
tion force depends on the density, which is different from many previous studies 
on kinetic-fluid models and is more physical in modeling but significantly more 
difficult in analysis. New approach and techniques are introduced to deal with 
the strong coupling of the fluid and the particles. The global well-posedness 
of strong solution in the three-dimensional whole space is established when 
the initial data is a small perturbation of some given equilibrium. Moreover, 
the algebraic rate of convergence of solution toward the equilibrium state is 
obtained. For the periodic domain the same global well-posedness result still 
holds while the convergence rate is exponential. 


1. Introduction 


1.1. The model. Kinetic-fluid models are widely used in the description of the 
dynamics of disperse two-phase flows. In these two-phase flows, the disperse phase 
is usually considered from the statistical point of view (kinetic equation) while the 
dense phase is from the hydrodynamic one (fluid equations). The kinetic equation 
is coupled with the fluid equations through the friction force. Typical applications 
of two-phase flows include the dynamics of sprays mm. diesel engines [5^15511551 
155 ] . pollution settling processes [ 5 ], rain formation [ 15 ], wastewater treatment [ 5 ], 
biomedical flows [2], combustion theory |36j . and so on. 

In this paper we are concerned with the following system of partial differential 
equations (see [TT]) of fluid-particle flows: 


dtF -b V ■ yxF = n\7v ■ [(u — u)F + VyF] 
dtn + V • (nu) = 0, 


( 1 . 1 ) 

( 1 . 2 ) 



(1.3) 
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with the initial data 


{F,n,u)\t=o = {Fo{x,v),no{x),uo{x)). 


(1.4) 


Here, the unknowns are F = F{t,x,v) > 0 for (t,x,v) € R"*" x fl x R^, denoting 
the density distribution function of particles in the phase space; and n = n(t, x) > 
0,u = u{t,x) G R^ for {t,x) G R"*" x fl, denoting the mass density and the velocity 
field respectively. The pressure function p depends only on n satisfying p'{-) > 0. 
In our present work, we take p{n) = Cqu'’' with the constants 7 > 1 and Cq > 0 
for simplicity. The spatial domain is H = R^ or (a periodic domain in R^). 
Compared with the model introduced in here we have normalized the phys¬ 
ical constants to be one for simplicity and added the viscous term —Ait in the 
momentum equation da. 

It is easy to check that for smooth solutions of the compressible Navier-Stokes- 
Vlasov-Fokker-Planck system (HD-dESl), the following quantities are conserved/ 
dissipated: 

• mass conservation: 


_d 

dt 



F dxdn = 0, 



0 , 


• momentum conservation: 

_d 
dt 


< / nudx+ I'T'dccdn > = 0, 

'■ Jn JJqxk^ ^ 


and energy/entropy dissipation: 


d 

dt, 


/ l^F^dnjdx-k J \Vu\‘^dx 


nxR3 


\{v-u)F-V^F\^ 

n - — -dxdi? 

F 


(1.5) 


with A = f’^ ^^^drj. 
Set 


M = M{v) = 


■ exp 




( 271 ) 3/2 

From the energy/entropy dissipation (11.511 . we know that {F, n, u) = (M, 1,0) is an 
equilibrium of the system (HUi-dESl). Thus it is natural to introduce the transforms 

F = M + y/Mf, n = l + p (1.6) 


to rewrite the system (HHl-dESl) as 

dtf + v- Va,/ -I- u ■ Vyf - i-u • w/ - M • v^/m 

= Ff + p(^Cf - u ■ V„/ + ^u-vf + u- vVm^ , (1.7) 

citp-I-M • Vp-I-(1-I-p)divM = 0, (1.8) 

dtu + u ■ Vu -f ——=- Au — uil -I- a) -I- 6. (1-9) 

1 -I- p 1 + P 
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Correspondingly, the initial data becomes 

/ Fq — M \ 

{f,n,u)\t=o = {fo{x,v),po{x),uo{x)) = ^ ,no(x) - l,uo(x)J. (1.10) 

In (|1.7I) - (I1.9I) . £ is the linearized Fokker-Planck operator defined by 


^f = 






and a = ,b = , depending on /, are the moments of / defined by 

a^{t,x)= / \/Mf{t,x,v) dv, b^{t,x)= / v\/M f{t, x,v) dv. 

Jr3 

1.2. Notations. Let ^{v) = 1 + and denote | • |j/ by 

\9\l-= [ + g = g{v). 

We use (•, •) to denote the inner product over the Hilbert space i.e., 

{g,h):= g{v)h{v)dv, g,h€Ll. 

JR^ 

For simplicity, we shall use || • || to denote the norm of or „ when there is no 
confusion. Define 

\\g\\l-=[[ [\yvg{x,v)\^+ iy{v)\g{x,v)\^]dxdv, g = g{x,v). 

For multi-indices a = ( 01 , 02 , 03 ) and P = (/3i,/I 2 ,/Sa), we denote by 

dp ^ 

the partial derivatives with respect to a; = (a;i,a: 2 ,X 3 ) and v = {vi,V 2 ,v^). The 
length of o and /3 are defined as |o| = 01 - 1 - 02 - 1-03 and \P\ = Pi + P 2 + Pa- We 
shall use the following norms: 

E ll^“5ll, \\g\\Hi^^ := E \\9pgl 

|C(|<S |a|-|-|/3|<s 

For g{t,x,v), we decompose it as the sum of the fluid part Pg and the particle 
part {I-Pjg: 

5 = P5 + {I-P}ff. (1.11) 

Here the projection operator P is defined as follows: 

P: —>• Span {a/m, ?;i\/m, 112 a/m, ?; 3 \/M}, 

and 

P:=Po©Pi, Pof: = a^VM, Pi/:= 6 ^ • uVM. 

This macro-micro decomposition is initiated by Guo [24] for the Boltzmann equation 
and later introduced in [16] to study the Fokker-Planck type equations. Notice that 
the operator C satisfies 

-[ g£gdv>\o\{l-Po}g\l, \Jg = g{v), 
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for some positive constant Ao > 0. Meanwhile, Cg can be computed as 
Cg = C{1 - P}g + CPg = £{I - P}g - P,g. 

Therefore, we have 

{-C{l-P}g,g)>Xo\{l-P}g\l, {-Cg,g) > \o\{l-P}g\l + m^. ( 1 . 12 ) 

For brevity, we define the temporal energy functional and the corresponding 
dissipation rate for (^f{t,x,v),p{t,x),u{t,x)) as the following: 


^o(i) := E E / - P}f)dx 


|o!|<3 i,j 


- E / dyd^y^-bdx, 

£i{t)-.=\\fr+\\pr+\\ur+ ^ ^iia“/f+ 

l<|a|<4 


(1.13) 


\/p'(l + p) 


1 + p 






+ Ti£o{t)+T2 E [ 9“M-9“Vpdx, 

|a|<3-^«^ 

:=l|V(a,&,p,u)||^3 + ||6-u||^4 

(1.14) 

+ E (||5“{I-P}/f. + l|5“Vuf), 

l“l<4 

(1.15) 

£2it):= E E IIW-P}/f, 

l<fc<4 \/3\=k 

|a| + |/3|<4 

(1.16) 

V2{t):= E IIW-P}/II- 

1<I/3|<4 
l“l + l/3|<4 

(1.17) 

£{t) :=£i{t) + T3£2{t), 

(1.18) 

'D{t) ■.= 'Di{t) + T3'D2{t), 

(1.19) 


where ri, r 2 , ra, C'fc(l < k < 4) are suitable constants to be chosen later. In addition, 
in torus we know the Poincare inequality is true, thus, the total dissipation rate is 
slightly different from 'D{t). We note that 

+ T4(||a|||2 + ||p|li 2 ) + T 5 II 6 + m||| 2 , 

where T 4 andrs are sufficiently small to be chosen later. Throughout this paper the 
letter C denotes a positive (generally large) constant and A a positive (generally 
small) canstant, where both C and A may change from line to line. The symbol 
A ^ B means < B < CA for some constant C > 0. 
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1.3. Main results. Our aim is to establish the global well-posedness and large¬ 
time behavior of strong solutions when the norm of the initial data ll/ollpf^^ + 
II (po) Wo)IIFT'i is sufficiently small (near the equilibrium). We also obtain the different 
time-decay rates depending on the spatial domain or . We now state the main 
results as follows. 

Theorem 1.1. Let 17 = and (fo, po,uo) be the initial data. Suppose that 
Fq = M + a/M/o > 0, and there exists a constant eo > 0 such that ||/o||ip4^ + 
IKpo, wo)||i/4 < eg- Then, the Cauchy problem p.7|l - (ll.l0ll admits a unique global 
solution {f,p,u) satisfying F = M + y/Mf > 0 and 

/ e C([0,oo);i7^(R3 X R^)); p, u G C{[0,oo); H‘^(R^)); 

sup(||/(t)||ff4 ^ -I- ||(p, u)(t)||// 4 ) < cdl/ollffi^ -I- ||(po,wo)||ff4), 
t>o 

for some constant C > 0. Moreover, if we further assume that 

+ WpmWm + l|w*(0)f^3 < +0O, (1.20) 

then 

I ll^|/lli= j + ^ + (1-21) 

^ |aH-|;S|<l 

for some constant C > 0 and all t > 0. 

Remark 1.1. In the assumption (11.201) . /t(0) is indeed defined through the Vlasov- 
Fokker-Planck equation (EB as follows: 

/t(0) ■=-v- Va;/o - Mo • V„/o -b ^mq • vfo + Uq ■ vy/ld 

+ Tfo + p(^T,fo - Mo ■ V^/o -b ^Mo • vfo + Mo • v\/m'^ ; 

and pt ( 0 ) and ut ( 0 ) are defined similarly. 

Theorem 1.2. Let 17 = and {fo, Pq,uo) be the initial data. Suppose that 
Fq = M + y/Mfo > 0, and there exists a constant eo > 0 such that ||/o||ii4 ^ + 
II(Po,mo)||h4 < eo, and 

/ aoda; = 0 , / poda; = 0 , / (6o + (1 + Po)mo) da; = 0 , 

Jt3 Jj 3 Jf3 

where ao = y/M foix,v) dv and bo = V3/M fo{x,v) dv. Then, the Cauchy 

problem (I1.7D - (I1.10D admits a unique global solution {f,p,u) satisfying F = M + 
'/Mf > 0 and 

f e C([ 0 ,oo);i 74 (T 3 X R^)); p,uG C{[d,^)-H\T^))- 

ll/(i)llil*„ + ll(P:M)(t)||_y4 < C{\\fo\\Ht^^ + ||po,Mo||ii'4)e“'^‘, 

for some constant A > 0 and any t > 0. 
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We remark that in the papers [T^ITBlfTT] for the related systems the time-decay 
rates are optimal in the whole space case. Here our time-decay rate (1 -I- in 
whole space case for the solution of (fTTD-ira is not optimal. The main reason is 
that, due to the strong coupling of the nonlinear terms in (ll.7l) - (ll.9l) . the spectral 
analysis cannot be carried out directly. We thus present another way to obtain 
the decay rate. In the periodic case, by the Poincare inequality we obtain the 
exponential decay which coincides with those in p^fTHlITT] in some sense. 

1.4. Some known results on kinetic-fluid models. There exist many versions 
or variants of kinetic-fluid models, depending on the physical regimes under con¬ 
sideration, such as the compressibility of the fluid, viscosity of the fluid, species of 
particles, interactions between the fluid and particles, motion of the particles, and 
so on. Below we review some kinetic-fluid models related to our system (fTB-lfOl). 
We discuss the case that both the fluid and the particle phases are isothermal. 
For the kinetic-fluid models with energy exchange involved, we refer the reader 
to [61E2]. Generally speaking, the kinetic-fluid models can be divided into two 
categories: incompressible models and compressible models. 

1.4.1. Incompressible kinetic-fluid models. If we assume that the fluid is incom¬ 
pressible, we obtain the incompressible kinetic-fluid models. The mathematical 
analysis of incompressible kinetic-fluid models has received much attention recently. 
In [27] , Hamdache established global existence and large-time behavior of solutions 
for the Vlasov-Stokes system. Boudin, Desvillettes, Grandmont and Moussa |7] 
proved the global existence of weak solutions to the incompressible Vlasov-Navier- 
Stokes system on a periodic domain. Later, this result was extended to a bounded 
domain by Yu m- Goudon, He, Moussa and Zhang m established the global 
existence of classical solutions near the equilibrium for the incompressible Navier- 
Stokes-Vlasov-Fokker-Planck system, meanwhile Garrillo, Duan and Moussa [12] 
studied the corresponding inviscid case. Ghae, Kang and Lee m obtained the 
global existence of weak and classical solutions for the Navier-Stokes-Vlasov-Fokker- 
Planck equations in a torus. Benjelloun, Desvillettes and Moussa [1] obtained 
the existence of global weak solutions to the incompressible Vlasov-Navier-Stokes 
system with a fragmentation kernel. Goudon, Jabin and Vasseur [201121] investi¬ 
gated the hydrodynamic limits to the incompressible Vlasov-Navier-Stokes system 
by means of some scaling and convergence methods. 

Assume that the fluid is incompressible and inhomogeneous, Wang and Yu [33] 
obtained the global weak solution to the Navier-Stokes-Vlasov equations, while 
Goudon, Jin, Liu and Yan [23] presented some numerical analysis on this model. 

1.4.2. Compressible kinetic-fluid models. Mellet and Vasseur [501 [31] studied the 
following compressible Navier-Stokes-Vlasov-Fokker-Planck system: 

dtf + v- V,/ + dW^iFdf - V„/) = 0, 
dtn -\- Vx • (nu) = 0, 

dt{nu) -f Va; • {nu 0 u) — Am -|- = — 



Fdfdv, 


( 1 . 22 ) 
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where 

Fd = Fo{u — v), Fo > 0 a constant. (1-23) 

In [30] they obtained the global existence of weak solutions of (I1.22I1 . and in |3T] 
they studied the asymptotic analysis of the solutions. In [14], Chae, Kang and Lee 
studied the existence of the global classical solutions close to an equilibrium and 
obtained exponential decay of the system to the system ( 11 . 221 ) . 

When the viscous term —Am in the system (11.221) is dropped, Duan and Liu [T^ 
studied the global well-posedenss of small solution in the perturbation framework. 
Carrillo and Goudon [10] investigated the dissipative quantities, equilibria and their 
stability properties. Morevoer, they also studied some asymptotic problems and the 
derivation of macroscopic two-phase models. 

As pointed out in [31] , the choice of drag force (11.231) may not be the most relevant 
one from a physical point of view. It could be more relevant from a physical point 
of view to assume that Fd depends on the density of the fluid, such as 

Fd = n{u-v). (1.24) 

To our best knowledge, the first rigorous mathematical result concerning the case of 
the drag force depending on the density n was obtained by Baranger and Desvillettes 
[3], where the following inviscid system; 

r dtf + v- V^f + nV„ • (/(m - v)) = 0, 

I dtn + Vj: ■ {nu) = 0 , ^ 5 ) 

dt{nu) + yX ■ {nu ^ u) + = —n / f{u — v) dn 

^ Js.3 

was considered and the local-in-time classical solutions were given. In the 
model (jl.l|l - ()1.3ll was introduced but no mathematical result was presented. In 
this paper, we study the global existence and large-time behavior of solutions to 
the model (|l.ll) - ()1.3|l . Our result is the first one on the global existence of the strong 
solution to the kinetic-fluid model when the drag force depends on the density. 

1.5. Strategy of the proofs of our main results. Compared with the studies on 
the compressible kinetic-fluid models in literature, the analysis of the system (HD- 
m or ([ni)-([r9i) is more complicated and difficult in mathematics, as explained 
below. 

Our approach is different from those used to obtain the global existence of clas¬ 
sical solutions in [T2l[l3l[T9] on the incompressible kinetic-fluid models and [BE] 
on the compressible kinetic-fluid models where the drag force is independent of the 
density. Due to the nonlinear terms caused by the strong coupling of / and p in the 
system ((ni)-([r9l), we cannot use the existing results on the Vlasov-Fokker-Planck 
system to obtain the regularity of /. Thus, we have to deal with the derivative of 
the particle velocity v in our arguments. 

Our main difficulties in obtaining the large time behavior of solutions come 
again from the strong coupling of / and p in the system (|1.7ll - (ll.9p . It prevents 
us from taking the advantage of the linearized spectral analysis to gain the rate 
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of convergence of solutions as in [EHTllIIllIIl]. To overcome these difficulties, we 
shall construct some novel functionals and adopt with modification some techniques 
in m to deal with the compressible Navier-Stokes equations and in for the 
Landau equation. Through the detailed analysis on the strong coupling terms of / 
and p we obtain the desired estimates. We believe that the methods developed in 
this paper can be applied to study the more complicated models in 0[22] , which is 
our forthcoming research project. 

The rest of paper is organized as follows. In Section 2, we shall establish the 
global existence of classical solutions to the problem (I1.7I) - (I1.10I) in the spatial do¬ 
main or T^. By the fine energy estimates we mainly use the local existence 

of strong solutions and continuum argument, motivated by [251126] for the Boltz¬ 
mann and Landau equations. In Section 3, by means of the energy estimates for 
the temporal derivative of the system and the Gronwall-type inequality, 

we eventually obtain the large time behavior. In Section 4, we shall establish the 
uniform a priori estimates with the aid of some energy functionals and correspond¬ 
ing dissipation rate. Although the uniform a priori estimates obtained in Section 4 
are needed in Sections 2 and 3, we shall present this lengthy part in the last section 
of the paper for the convenience of readers. 

In the rest of this paper, we shall omit the integral domain x or in the 
integrals for simplicity. 


2. Global existence of the classical solutions 


In this section, we shall establish the global existence of classical solutions to 
the problem (ll.7ll - (|1.10ll in the spatial domain or T^. It is well known that 

by the uniform a priori estimates we shall obtain the global existence of solutions 
with the help of the local existence as well as the continuum argument, under the 
smallness and regularity conditions on the initial data. Here we first construct the 
iteration process to obtain the unique local solution, then the global existence of 
solutions follows from the continuum argument and the uniform a priori estimates 
obtained in Section 4. 

Now we define iteratively the sequence {F^, as the solutions to the 

system: 


' dtF^+^ -H V ■ V^F^+^ -{1 + -P V„F"+i) 


= -(I + p"K-V„F"+\ 

^ -f u" • -f (1 -f p^) div u” = 0, 



Setting F^ = M + we can rewrite the above system as 

+V ■ v,r+i - (1+p")£r+i 



( 2 . 1 ) 
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+ w" • Vp’^+1 + (1 + p") divM”+l = 0 
1 


( 2 . 2 ) 


9tu"+i-Au”+i = -u” • Vu" + 7(1 + p^y-'^Vp^ + u"(l + a") + 6", 

1 + p" 

(2.3) 

where n = 0,1, 2 ,..., and (/°, p°, u°) = (/o, Po,uq) is the starting value of iteration. 
We define the solution space X{0,T-,A) by 

f G C{[0,T],H'^{nxM.y), M+ yMf >0; 
p G c([o, TiHyn)) n ci([o, TiHyn))- 
u G C([0, T],H\n)) n Cl([0, T],H\n))- 
supo<t<T{||/Wllif^.„ + II(p,u)||h4} < A; 
pi = ^(-1 +infp(0,a;)) > -1; 

info<t<T.a:eop(i,ai) > pi. 

The main result of this section reads as follows. 


XiO,T;A) := 


(2.4) 


Theorem 2.1. There exist Aq > 0 and T* > 0, sueh that if fo G x ]R^),po G 

H‘^{Q),uo G with Fq = M + a/M/o > 0 and £{0) < with £’(0) ^ 

+ IKpoi wo)||^4 (see the specific definition of £{0) in Section 4), then for 
each n > 1, {f'^,p^,uy is well-defined with 

(r,p",n")G A(0 ,T*;AIo). (2.5) 

Moreover, the following statements hold: 

(1) (/^p^w") n>o is a Cauchy sequence in the Banach space C([0, T*]; iJ^(n x 
IR3)) X c([o,r*],ij3(fi)) X c([o,r*],ij3(f2)), 

(2) the corresponding limit function denoted by {f,p,u) belong to X{0,T*] Aq), 

(3) {f,p,u) are solutions to the Cauchy problem (I1.7I) - (I1.10|) . 

(4) {f,p,u) is unique in X(Q,T*-, Afi) for the problem (I1.7I) - (I1.10I) . 

Proof. Let T* > 0 be a constant which will be fixed later. For brevity, we can 
assume that {f^,p^,uy are smooth enough in order to take the forthcoming cal¬ 
culations, otherwise, we can consider the following regularized iteration system: 

^ - (1 -b p”’'^)£/”+i’'= 

= -(1 + p^yu^^^ ■ - v^/m) , 

dtp^+^’^ + • Vp"+i''= -b (1 -b p"’'^) div = 0, 

- 1 - Au^+^'^ = -U^'^ ■ 

1 -b p"’'" 

+ 7(1 + p^'y-'^xp^'^ + u’^’yi + ay + 

(/”+'’% «”+'’')(0) = (/o%Po,«o) 

for any e > 0 with (/o,Po,Mo) a smooth approximation of {fo,po,uo) and pass to 
the limit by letting e —>■ 0. 
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Applying 9“ with |a| < 4 to the equation m, multiplying the result by 9“/"+^ 
and then taking integration over one has 

1 ^pa/n+i||2 + J{i + p")(-/:9“/"+i,a“r+i) dx 

= / d^p^{cd°‘-^r+\d‘^r+^)dx 

a“{(l + + yM/”+^)}5“/”+^ da;di; 

<C(l + ||p'^||^4)|K||^4||r+i||^.(^4) 

+c(i + iip"ii^4)ii«-|i^4||/"+i^.(^4)iia“r+i. 

+ c\\p-\\h^{ ^ (2.6) 

|a'|<3 

Notice that 

J{i+ p")(-£a“r+i, a“r+i) dx > A|i{i - Po}d^r+^\\i. 

By adding \\Pod °‘to both sides on the inequality (12.6|) . then, taking sum¬ 
mation over |a| < 4, we have 

ii5“r+i^ + A ^ \\d^r+x 

\a\<4 

<C(l+||p"||^4)||ll-||^4||r+l|l^.(^4) 

+ C'(l + ||p"||?,4)||«”||?,4||r+l|li.(^4) 

+ C\\p-\\h- E WdV^^X + CWr^TLliH^y 

|a|<4 

Similarly, for any 0 <t <T < T*, we obtain that 

|a| + |/3|<4 

< c(i + +cwp^h^ e 

|a| + |/3|<4 

+ ^(1 + \\p-\\u\w^ma + wr^Tmj- ( 2 - 7 ) 

Next, according to [29], for the system (12.2p and ()2.3|) . there exists a unique 
solution satisfying p"+^ > pi, and 

p^+^ G ci[o,T],H\n))nc\[o,T],H^n)), 

y"+i g c{[o,T],H‘^(n))nc\[o,T],H^{n)). 

Now we estimate ^"^^)lll 44 . Applying d°‘ (|a| < 4) to the system (12.21) 

and (1^ . multiplying the results by and respectively, and then 

taking integration and summation, one has 



1 d ^ 
2dt ^ 


|a|<4 
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< (7(1 + ||p"'|||/4 + ||'u”|||^4)||p"^^|||f-4 + (7(1 + ||p"||fl-4)||'«"''"^|||f-4 

+ C'(h"||^4 + Iir||?,4)(l + \\u-\\l.) + C\\p-\\\.{1 + \\p-\\%.). (2.8) 

Adding up (I2.7p and ()2.8p gives 

+ A ^ lia^r+'f. + A ^ ||V5“«-+if 

l“l + li®l<4 |a|<4 

< c(i + iip"iii,4)iir+ii,4^^+ c{i + iip'^ii^4)ii«'^iii,4(i + wr^Ymj 

+ C{1 + ||p"f^4 + \Krm)\\p-+Ym + C{1 + ||p"f^4)||n"+i?,4 

+ ^7(1 + ||«"||l,4)(||rf^4^^ + ||n"f^4) + C\\p-\\U^ + \\p-\\%.) 

+ C\\p-\\h^ J2 mr+X- ( 2 . 9 ) 

|a| + |/3|<4 

Using induction, we may assume A„(T) < Aq and A„(0) < X some Aq > 0 
with 

An{T):= sup {|| p ”( t )||^4 + || M ”( t )|| 1^4 + ||/”( t )||^4 }. 

0<t<T ’ 

Integrating (1^ over [0,r] gives 

A„+i(r) + A r{ J2 \\yd^u-+X+ J2 \\dX^X\l}dt 

l«l<4 |a| + |/3|<4 

< A„+i(0) + C(1 + a| (T) + A2 (r))A„+i(T)r + C(A„(r) + Ai{T))T 

+ cAk{T) J2 

|a| + |/3|<4 

< ^ + C(1 + Al)TAn+i{T) + C{Ao + Al)T 

+ CAl{T) Y. \\dprX\l- ( 2 . 10 ) 

|a| + |/3|<4 

It follows that, for T <T*, 

(l — (7(1 + Aq)T)A„+i(T) < ——h C{Ao + Aq)T. 

Choosing T* satisfying T* < ^, and Aq sufficiently small, we conclude that 

An+l < Aq. 

For the equation of with the help of the maximum principle, we have 

= M + > 0. 

Now we explain that ||/"^^||^4 is continuous over 0 < t < T*. In fact, it follows 
from the inequality: 
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= 1 

<CaI Y. t\\d^r^Y.dv + C{Ao + Al)\t-s\, (2.11) 

|a| + |/3|<4-^'* 

which can be proved by the same process as the proof of (IQ) . Meanwhile, 
||^ayn+i ||2 jg j^tegrable over [0,T*]. Hence, (j2.5ll holds true for n + 1 and so 
it does for any n > 0. 

Next, we study the following system: 

-n+v Vx(r+' - n - - n 

= p^cr+^ + (1 + - V„r+i) 

- (1+ p^-^)u'^-^{vVm+- v„r) - p^-^cr, 

dt(u^+^ - It")-A(w"+1 - u") 

^ ' 1 + p^ ^ ’ 

= (t^ - 1 

Vl + p" l + p"-i/ 

- ((tt" - u"-i) • Vit" + • V(it" - u"-i)) 

+ 7((1 + p")^-"Vp" - (1 + p"-i)'>'-2Vp"-^) 

+ 6" _ + (a'l _ a"-i)rt"-i + (n" - rt"-^)(l + a"), 

ai(p"+i - p") + u" • V(p"+i - p") 

= -(u" - u"-i) • Vp" - (1 + p") du(tt"+i - u") - (p" - p"-i) divrt". 
Similarly to (I2.9L we obtain that 

- n\Hi^ + - P^Wm + 

|a| + |P|<3 

<C(l+||p"|||,. + ||zt"r^.) 

X (11/"+' - .r||?^a+ ||p"+' - P"||l,3 + |k"+' - n"||^3) 

+ C||p"||^3 ||a“(/"+l-/")||2+C(l + |K-'||?,3)||/"-/"-'|l?,3 

|a| + |P|<3 

+ C^{(1 + I|p”ll?t4)(l + ||/"||1,3) + ||(u"-',/",u")||l,4}||u" - U"-'||?,3 

+ c|l + ||(n",n"-')||l,4 + Up\u-)\\U\plfH^ + (1 + h”-'ll?t3)||/"|H,3 

+ E \\dpn\iW-p^-Tm- 

l«l + l/3|<3 ^ 
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Here we have used the Sobolev embedding ^ Since f (0), T*, and 

Aq are sufficiently small, by using (I2.10L we know that 



is also sufficiently small. Hence, there exists a constant k < 1, such that 
sup {iir+i - + iiu"+i - 

0<t<T* 

<K sup {||/"-/-i ||^4 + ||p"-p"-i ||^4 + ||n"-n"-i||^ 4 }. ( 2 . 12 ) 

0<t<T* 

According to p.l2(l . we conclude that (/",p",'u”)„>o is a Cauchy sequence in 
the Banach space C{[0,T*], x R^)) x C([0,T*],i73(H)) x C{[0,T%H^{n)). 
Hence, in this Banach space, there exists a limit function (/, p, u) such that (/, p, u) 
is a solution to the Cauchy problem (I1.7I) - (I1.10I) by letting n —oo. From the fact 
that F^{t,x,v) > 0 and the Sobolev embedding theorem, we deduce that 

F{t,x,v)>0, sup ||/(t )||^4 <Ao. 

0<t<T* 

Similarly to the proof of (j2.11ll . we see that / £ C[[0,T*], H^{n x R^)). Thus, we 
can conclude that {f,p,u) £ X{0,T*; Aq). Finally, let [J,p,u) £ A(0,T*,Ao) be 
another solution to the Cauchy problem (I1.7I) - (I1.10I) . By taking the similar process 
to that of ( 12 . 121 ) . we have 

sup { 11 / - /lli/4 + Up - p||^4 + ||m - u \\ h ^] 

0<t<T* 

<K sup (11/-/||ii4 + Up - p||j^4 + ||ii- 1211 ^^ 4 } 

0<i<T* 

for K < 1. Hence we deduce that f = f, p = p^u = u and uniqueness follows. □ 

Since £1(0) ||/o||ff 4 ^ + ||(po,wo)||ff 4 , there exists cq > 0, such that if ||/o||ii 4 „ + 

IKpo, uo )||//4 < Co, we have £1(0) < Next, with the aid of Theorem 12.11 we 
obtain the local solution on [0, T*] that satisfies the uniform a priori estimate (j4.28l) 
in Section 4. Finally, by taking the standard bootstrap arguments similar to those 
in [T6l[24l[^, we obtain the global existence and uniqueness of classical solutions 
in both Theorem 11.11 and Theorem 11.21 


3. Large time behavior of the classical solutions 

In this section, we investigate the time-decay rates of global solutions to the 
problem (I1.7I) - (I1.10I) . We obtain the algebraic rate of convergence of solution toward 
the equilibrium state in the whole space case, while for the periodic domain case, 
the convergence rate is exponential. 
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3.1. The case of the whole space. In this subsection we consider the large time 
behavior of classical solutions in the whole space 17 = In order to obtain the 
desired decay rate (11.211) in Theorem ll.il we first introduce some new functionals 
which are similar to those |151I28) in spirit. Then we perform the energy analysis 
to the temporal derivative of the system (inii-diii), instead of the original equa¬ 
tions, to gain the one-order derivative. Finally, we combine the energy estimates 
together with the uniform a priori estimates obtained in Section 4 and the following 
Gronwall-type inequality to obtain the time decay of the solutions. For simplicity, 
we shall denote {ft,pt,ut) = {dtf,dtp,dtu). 


Lemma 3.1 ( [H]). Let y{t) G C'^([to,oo)) satisfy y{t) > 0,T = /j“?/(s)ds < -l-oo 
and y'{t) < a{t)y{t) for all t > to. If a{t) > 0 and B = a(s)ds < -l-oo, then 


< [hyitp) + 1) exp(A + B)-l 


First of all, we consider the system for {ft,pt,ut). 


dtft + V ■ V^ft - (1 + p)Cft = PtCf - (1 + p)u ■ {Vvft - 

- [ptu -h (1 -H p)ut) ■ (Vvf - , 

dtpt + u ■ Vpt + pt divM = -UtVp - (1 -I- p) divMj, 

dtut + Ufyu + u- S/ut - Aut -I- ^ Vpt -I- ^ 

1 -I- p 1 + P (1 + P) 

= - 7(7 - 2)(1 -I- py~^ptVp - ut{l -I- a) - uat + bt. 


Pt An 


(3.1) 

(3.2) 


(3.3) 


Notice that there are some similar structures between the system dSU-dSSl) and 
the system (ll.7l) - (ll.9l) . and we can benefit from the proofs of Lemmas l4.1ll4A] Thus, 
we will omit some details of the proof of the following lemma. 


Proposition 3.1. Assume that {f,p,u) is the solution obtained in Theorem \1.1[ 
Then, we have 

2 + llPtiP + ~ ~ ^tll^ 

< eWVipt, btW + C.(||{I - P}f\\l + ||u - 6f + II Vull^O 
+ ^(l + ||p||^.)||u||^.||{I - + C{1 + IIpII?,. + ||uf^.) 

^ l«l<2 ^ 

X iwftr + upuutnj,.) (3.4) 

for e > 0 sufficiently small. 


Proof. First, multiplying (13.1() by ft, (|3.2|) by pt, and ()3.3I1 by ut, respectively, and 
then integrating over R.^ and taking summation, we obtain that 

2^(ll/ill^ + llPtll^ + IktlP) -I- A||{I - P}ft\\l + AllVutlP + ||ut - ^fll^ 
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< JJ ptCfftdxdv + j ptu-btdx- JJ {ptu + {1 + p)ut) ■ (v v - f ft dxdv 

— J au^ dx + — JJ (1 + p)u ■ vf^ dxdi; — J atu ■ utdx + J p{ut — h) ■ bt dx 

- J Uf Vppt dx + 7(7 - 2) y (1 + pV~^PtUt ■ Vpdx 

- J {I + p) div utpt ~ J ^ 

— J u - y ptPt ~ J Pt udx — J Ut ■ Vu ■ Utdx — J u- \7ut ■ Ut dx 


1 


(1 + p)' 


-Vp • Vmj ■ Utdx — 


1 


-ptAw • Ut dx. 


(3.5) 


(1 + p)^' 

Next, we estimate each term on the right hand side of (I3.5|) . For the first two terms, 
we rewrite them as follows: 


PtCf ftdxdv + J 

ptn ■ bt da 



-DC- 

-P}/(v 


1 

+ 

■p}/(^ 



2 





ptu ■ bt dx 


Hi + 112 + Ila + 114. 
For the term IIi, one has 


(3.6) 


Hi <c\\pt\\ J Y. + + 


l«l <2 


LI 


dx 


<ciiptii E ii5“{i-p}/iuii{i-p}/tiu 

l«l <2 

<e\\{l-P}ft\\l + CM^ E l|5“{I-P}/ll^. 

I«l <2 

By the definition of P and noticing the equality (V^ + f )ff = b^'/M for any g, one 
has 


JJ Pt {I - P}/ • btVM dxdx 


n 2 

<C||V&*||^i||p*||||{I-P}/||. 
<e||V5*||U+C',||{I-P}/||^|lp*f 
for any e > 0. Similarly, for the terms Ila and 114 , we have 


n. 


JJ Pt (V„ + ^) {I - P}/ • bt^/Mdxdv 
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<C||V 6 t||^i||pt||||{I-P}/||. 

<e\\{i-p}M\l + c,\\m\\U\ptr, 

114 = J pt{u-b) ■ btdx 
<C\\Vbt\\Hi\\u-b\\\\p,\\ 

<e\\vh\\],,+c 4 u-b\mp,r 

for any e > 0. Plugging the above estimates into (13.61) . we can rewrite it as 

JJ pt£fftdxdv + J ptu-btdx<e{\\Vbt\\l^ + \\{l-P}ft\\l) 
+c,( ^ \\d‘^{i-p}f\\i+\\u-br+\\vbrm)\\ptr 

\ a \<2 

for any e > 0 . 

For the third and fourth terms on the right hand side of (13.51) . we have 
+ (1 + p)ut) ■ fft dxdv - J auf da; 

, + ^yi-P}ftdxdv 


(3.7) 


V„ + 


|)p/. 


dxdu — / auf da;. 


(3.8) 


We have the following estimates for the right hand side terms in (13.81) : 


„ + - P}/tda;du 

<cj iWufh^WptW + ||/IU~(i + llp||L<»)ht||)||(v. + |){i-p}/,|idu 

< C (1 + WpWh^ + h||4a2)||(pt,«0lll|Vx/|U2(^x)||{I - P}M\. 

<e\\{l-P}M\l + C,{l + \\{p,u)rH.) 

X ( ||9“{I-P}/f + ^ ||V4a,6)f)||(p,,u0f, 

l«l<2 |a|<l 

JJ {ptu + [1 + p)ut) f (y V + JjPft dxdv - J auf dx 
= J {ptu + put)abt + J aut ■ (bt - ut) dx 

< e\\ut - btr + a||Vaf^.(||/,|p + ||u*||2) + C.IKVp, Vu)||^.(||p*|p + ) 

for any e > 0. Plugging the above estimates into (13.81) gives 

- JJ {ptu + (1 + p)ut) (yy-Jjf-ft da;du - J auf dx 

< e(||Mt - &t|p + ||{I - P}ft\\l) + Ce{l + \\p\\']j2 + IIMIII42) 

{ ^ ||a“{I-P}/f+ ||V(a,&,p,u)||^,}(|l/*f+ ||(pt,u*)f)- (3.9) 

l«l<2 


X 
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The remaining terms on right hand side of ()3.5I1 can be treated easily. Hence we 
only list the bounds below and omit the details for brevity. We have 

JJ 2 P')^ ' daidi; — J atu ■ Ut da; 

< <\\ut - hf + ||{I - P}/*f) + C,{(1 + ||p||l,.)ll Vu||^. 

+ ||Vp||^i||Vir||^i}||/tf+ C(l + ||p||ff.)||ir||^.||{I-P}A||2, (3.10) 

J p{ut - bt) • fotda; < e||Mt - + <^£11 Vp||^i ||/t|p, (3.11) 

y (1 + p) diYUtPt + y 

<4Vutf + C,\\Vprm\\pt\\\ 

y u-VptPt dx + J Pt divwda; < e||Vpt|p + CellVull^illptlP, 

J PtUt • Vpda: - 7(7 - 2) J {1 + p)'^~^ptUt ■ \ 7 pdx 
<e\\yutf + c4yp\\U\ptf, 

Jufyu-utdx + J u-Vut ■ utdx < e\\Vut\\'^ + Ce\\yu\\'jji\\ut\\'^, (3.15) 


j\7p ■ \7ut ■ utdx — 


jPtAut ■ Ut dx 


{l + P? - ' ' J (l+p)2' 

<e\\Vutf + C4{Vp,Vu)\\l4{pt.UtW. 


(3.12) 

(3.13) 

(3.14) 


(3.16) 


Plugging all the above estimates (13.71) . (13.91) . and p.lOI) - p.l6l) into (13.51) and then 
choosing e sufficiently small yield (13.41) . □ 


Proposition 3.2. Assume that {f,p,u) is the solution obtained in Theorem M. 1 
Then, we have 




l<|a|<3 


Vp'i^ + p) 2 ^ 


1+P 


+\ ^ {\\d-{i-p}f4i+\\d4ut-btW + \\vd‘^utr} 

l<\a\<3 

< e(||Vat||^i + ||V&t||^2) +C^{l + ||p||//3)||Vu||^2||/t||^2(^3) 

+ cJ ll^^“{I-P}/f. + l|V.(a,5,p,u)||^2|||(pt,Ut)||?,3 

^ |a|<3 ^ 

+ C'l|p||ff3||V(at, bt)\\‘jj2 + (7(1 + ||p||Ef3)||u||^4||Vpt||^2 
+ C{(l + ||p||^3)||u||^3+||p||^3} ^ \\d‘^{I-P}ft\\l 

l<|a|<3 

+ c| ^ - P}/f. + II V(a, b, p)||^2 + (1 + ||p||?^3)|| Vulll,. |>||/t|li3(773) 

^ l«l<3 
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+ C' E ||a“{I-P}/f.+ l|V(a,6,p,«)||^3 + ||Vwf^4 

^ |a|<3 

for e > 0 sufficiently small. 


(3.17) 


Proof. Applying 9“(1 < |a| < 3) to (I3.1D - (I3.3I) . multiplying the results by d°‘ft, 
\iXp)^ and 9“ut respectively, and then adding them together, we obtain 


+ J{-c{i-p}d^ft,d^ft)dx + \\d^{ut-bt)\\^ + y^|va“ut|2dx 

: JJ 8“ (^{ptu + (1 + p)ut)^f^d°ftdxdv - JJ d°{{ptu + (1 + p)ut) ■ V„/)5“/t dxdu 
+ JJ 5“ ((1 + p)u ■ ^ft) d^ft dxdv - j d^{uat) ■ d^^ut dx 

- y (1+p)'“ ■ + y da; 

- 1 d‘^{pcft)d-ftdx- 1 d^ 

+ y d“[ptU +put) ■ d'^btdx + j d°'ut ■ d'^pt dx 

+ y[-5“,pV2,>tLli±^a“pt da; + J [-8°^, u ■ d°‘pt dx 


1 f a ( P'i^ + P) 

(1 + p)^ 

oa 


(i + p)^ 

p'i^ + p). 


)\8°'ptf dx - J divu + Vp- 8°‘pt dx 


- u- SJd^ut ■ 8°‘ut dx- V 


1 + p 


(1 + p)^ 

V9“ut • 9“ tit dx - y 5“(ut • Vu) • a“ui dx 
- E 1 8^{j^)8^-'^Auf8^utdx- J 8^{-^j^^PtAu)-8^utdx 

pt ■ 8°‘ut dx 


1<7<q; 

+ [ [-8°‘,u-V]uf8°‘utdx + 


_ jAildv 

1 + p 


- 7(7 - 2) y a“((l + py-^ptVp) ■ 8^ut dx- J 8^{aut) ■ 8^ut dx. (3.18) 

Now we estiamte the terms on the right hand side of (13.181) . First, we have 

JJ 8^{iptU + put)-VJ)d^ftdxdv + J a“((ptu + pu0-^/)a“/idxdx 


< 


< 


, + 2 ) ^ 


HI 


Wd^fthldv 


C{\\yu\\H4Pt\\H^ + l|Vp||ff3||ut|U3)||(v„ + ^)Vx/y^ ||a“/t|| 
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<C{ ^ ||5“{I-P}/||2 + ||V(a,6)||^.}||9“/if+ q|(Vp,V«)||2,.||(p*,,xO||l,3, 

|a|<3 

JJ d°‘{ut(yy + ^^f)d“ftdxdv 

<C||VMt||ff2||(v, + ^)v,/ 


L2(i/2) 


l|5“/t|| 


77||V,Xt||^.+C,{ ^ ||9“{I-P}/||^ + ||V(a,&)||?,.}||9“/t|p 

|a|<3 


with rj small enough. Thus the first two terms on the right hand side can be 
bounded by 

{ ^ ||9“{I-P}/||2 + ||V(a,6)||^.}||a“/tf+ C||(Vp,VM)||l,.||(pt,^rO||?^3 

|a|<3 

+ n\\Vut\\l.+C^{ ^ ||9“{I-P}/||2 + ||V(a,6)||^.}||a“/tf. (3.19) 

|a|<3 

For the third and forth terms on the right hand side of (I3.18|) . we can rewrite 
them as 


JJ d‘^(^il + p)u-^ft)d‘^ftdxdv- J d'^{uat)-d‘^utdx 

= JJ{1 + p)u ■ ^d'^ftd^ft dxd^; - J d^atu ■ d^ut dx 

+ E + p)u) ■ ^d^-ytd^ftdxdv- j d^-"<atd^u-d^utdx) 

= i y(l+p)u(?;, |{I-P}i9“/tp)dx + J d^atu- d°'ibt- ut)dx 

+ y’(l + p)u-(^;Pa“/t,{I-P}a“/t)dx + J pd^atu-d^bt dx 

+ E C^,^(^JJ d''{{l + p)u)-^d‘^-ytd‘^ftdxdv-J d^-^atd^u-d^utdxy 

1<7<q; 


For the terms on the right hand side of the above equality, we have 

ijil + p)u- (u, |{I - P}a“/tp) dx < C(1 + ||Vp||^3)||Vu||^3||{I - P}5“/t||2, 

J {1 + p)u ■ (xP9“/t, {I - P}a“/t) dx < C||(l + p)u\\l<^ ||9“/t||||{I - P}a“/t|| 

< ,7||{i - p}9“/*f + c,(i + ||p||l,3)|| vuiiu ||a“Af, 

J ud°'at ■ d°‘{bt - ut)dx < CH V-«||i/i ||5“at|| ||9“(ut - bt)\\ 

<v\\d‘^{b,-u,)r + cjvu\\Uft\\lump 
J pd°‘atu- d'^btdx < C\\pu\\L--\\d°‘at\\\\d°‘bt\\ < C(||Vp||U + \\'^u\\m)\\ft\\li(Hi)^ 
J d'^ud°'~'*at ■ d°‘ut dx < C'||Vu||ij 2 ||Vat||iji ||Mt||/i-3 
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< e||Vatlllfi + CellVu|||t-2llwt||lf-3, 

Jj d^{{l + p)u) ■ ^9“-'^/t9“/tdxd« 

= J d-^{{i + p)u)-^(d^-^i-P}ft + d^-^Pft)d^ftdxdv 

<c'(i + ||Hk3)||Va||^3{ ^ ||a“{i-p}/t||. + ||v(at,Mlk3}||a“/t|| 

l«l<2 

<77 ^ ||a“{I-P}/t||. + C,(l + ||p||^3)||Vaf^2|l/t||i3(^3) 
l«l<2 

+ e|| V(at, &t) 11^1 + (^^(l + ||p||^3)|| Vm||^2 W/tW^KH^)^ 
where e, 77 are sufficiently small constant. Thus the third and forth terms have the 
following bound: 

JJ d'^{{l+p)u ft dxdv- J d‘^iuat)-d‘^utdx 

< C{1 + \\yp\\H.)\\Vu\\H.\\{l-P}d‘^ft\\l + C||V(p,a)f^t||/t||i3(^3) 

+ V{J2 - ^‘)ll'} + G(1 + Ilpllff3)|| Vu||^2||/t||i3(^3) 

| a |<2 

+ e||V(at,6t)||l,t+C,(l + ||p||l,3)||Vu||^2||/t||i3(^3)+C',||Vu||^2||ut||?,3. (3.20) 

Similarly, we have 

(1 + P)'^ ■ V„]/t9“/t dxd-y 

<77 ^ ||a“{I-P}/t||. + C,(l + ||p||l,3)||Va||l,2||/t||i3(ff3) 
fol<2 

+ e||V(at, 6t)||^i +Ce{l + ||p||f73)||VM||^2|l/t||i2(jt'3). (3.21) 

Next, we deal with JJ d°'{pt£f)d°‘ftdxdv. We have 


JJ d°^{ptCf)d^ftdxdv 

= -JJ d‘^(^Pt{W^ + ^)f) -d^^btVMdxdv 

- II 5“(pt(V„ + ^)/) • (v„ + ^)5“{I-P}/7da:d7; 

< C\\Wpt\\H^ I (v„ + (||^^“{I - P}/tll. + ||5“fot||) 


't|lff3 


77||5“{I-P}/t||2+C,{ ^ ||a“{I-P}/ll" + l|V(a,l 

|a|<3 

+ e||V 6 tf^ 2 +C,{ ^ ||5“{I-P}/f, + ||V(a,6)||^2}||pt|ll,3. (3.22) 


|a|<3 


Similarly, we have 

- [[ d°‘{pCft)d^ftdxdv 
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= 11 + ^)ft) ■ (v. + ^)d:ftdxdv 

<C\\Vp\\^.{ Y, ||a“{I-P}/t||2 +||V(at,6011^2}. (3.23) 

|a|<3 


For the sake of brevity, we only give the bound of the remaining terms on the 
right hand side of (I3.18|) as follows: 


J d°‘(ptu + put) ■ d°‘btdx < e\\Vbt\\H2 + C'e||V(p, u)||^2|l(pt, Ut)||^ 3 , 

J 1 + P J K^ + P) 

<v\\yut\\l2 +c,\\wp\\l2\\pt\\l^. 

<C\\Vp\\H4^u\\H4dytf + C\\u\\H4ptfm, 

J[-d^,u-V]ptj^^^^d^ptdx- J d‘^{pt dYvu + UfVp) ^^^^J^ 

< »7||VMt|||^2 + C't7||Vp|||^3||Vpt||lf-2 + C'||M||/i-4||Vpt||l/2, 

. 1 


(3.24) 

(3.25) 

(3.26) 
d°‘pt dx 

(3.27) 


u ■ Vd°‘ut ■ d°'ut dx + 


1 + p 


Vd^^ut ■ d‘^ut dx 


+ Y Co.,1 j d^(-^)^°'~^^'^td°^'^tdx 

l< 7 <a ^ 

< + C^(||V(p,^)||?,. + ||Vpf^.)IKII?/3, 

J -d°‘{ut ■ Vw)a“ dx + y [-5“, w • V]utd°‘ut dx 

J L 1+P J 

< p||Vw*||^2 +C'^(||(Vp,Vu)||?,3 + \\VpfH4\\iPuUt)fH^. 
|9“ptp dx < ^(1 + ||p||^3)||«||^3|| Vp*||^ 2 , 


J ^“( (i+\) 2 PAn) -a^ntdx 


(1 + p)^' 

< Oi + llp|ll^3)||Vp||^2||pt||?,3 + c\\Vu\\U\u4m, 

a“((i + p)'^"^PtVp) • a“utdx < c(i + ||p||ff3)||Vp||^3||(pt,Mi)llff3, 
d°'{aut) ■ d°‘ut dx < p||VMt ||^2 + C'^||Va||^2||Mt||^3. 


(3.28) 


(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 


Finally, plugging the above estimates (I3.19I) - (I3.33I) into (13.181) . using (ll.l2l) . and 
choosing p sufficiently small, we obtain (13.171) . □ 
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Remark 3.1. Combining ProDOsitioii l3.ll and ProDOsition l3.21 for the solution 
{f,p,u) obtained in Theorem ll.il we have 


l_d 

2 dt 


|a|<3 |a|<3 l<|a|<3 




+ A E {||5“{I-P}/t|l2 + ||va“ixtf+ ||5“K-Mf} 

|a|<3 

<e{||Vat||^i + ||V 6 t ||^2 + ||Vpt|p} + Ce(l + ||p||/i- 3 )||Vu||^ 2 ||/t||i 2 (j/ 3 ) 

+ cJ E ll^^“{I-P}/ll" + l|V.(a,6,p,ix)||^2 + ||«-5f|||(p,,ii,)||^3 

^ |a|<3 ^ 

+ C'l|p||ff3||V(at, bt)\\'jj2 + (7(1 + ||p||ff3)||u||^f4||Vpt||^2 

+ C((l + ||p||^ 2 )h||^ 2 +||p||^ 3 ) E 

l<\a\<3 

+ C(l + ||p||^4 + ||u||^ 2 )| E l|a“{I-P}/||^+||V(a,&)||l ,2 

^ |a|<3 

+ Il^p|lff3 + l|V'u||^4|(||/t||^2(^3) + ||Pt,'Wt||^3). (3.34) 


Next, we need to estimate ||V(a 4 From (14.1211 - (14.141) . we deduce that 

{at, bt) satisfy the following system 


Stot + div6t = 0, (3.35) 

dtbt,i + diat + Ea,r„({i-p}/t) 

3 

= “(1 + p)bt,i — Ptbi + ptUi{l + a) + (1 + p)ut^i{l + a) + (1 + p)uiat, (3.36) 
^ibt,j 4“ ^jbtji (f 4” P^{^t,ibj 4“ Utpbi T utbtp Ujbt^i^ Pt{^ibj 4“ Ujbt^ 

= -9tr({I - P}ft) + r,., {It + rt+ St), (3.37) 


where 


It := -V ■ V,{I - P}ft + L{1 - P}ft, 

Tt :=at(-u-V4I-P}/+iu-u{I-P}/), 

St :=at(pM-^V„- (^VM{I-P}/ + ^{I-P}/-uyM{I-P}/)). 
Define the following functional 





({I-P}/t)dx 


(3.38) 


Similarly to the proofs of ProDOsition l4.3l and ProDOsition l4.41 we have the following 
two results: 
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Proposition 3.3. Assume that {f,p,u) is the solution obtained in Theorem \l. 1\ 
Then, we have 

<C[\\ut - bt\\‘\j 2 + ||{I - P+ C'(l + WpW’h^ + 11 ^ 11 ^ 3 ) 

X (II V(a, b, u)||^2 + ||{I — P}/||i2(//3) + ||w — &II// 2 ) {W/tWlKH^) + ll(Po wt)||^3) 
+ C'(||(p, m)||^ 3 + ||p||ff3||M||^3) {\\ut — bt\\'jj2 + ||{I — P}/t||^2(^3)). (3.39) 

Proposition 3.4. Assume that {f,p,u) is the solution obtained in Theorem \l.l\ 
Then, we have 

^ J d^utd^wptdx + ^p'{i)\\vd^pt\\^ 

< C(||Vut||^3 + ||ut - ^* 11 ^ 2 ) + C(l + 11 ^ 11 ^ 3 ) 

(ll^(OjP)llff2 + ll^^llu^) (||/t||i2(^3) + ||(Pt,Mt)||^3) 

+ C||(p,u)||^3||V(p*,Mt)||^2. (3.40) 


Now, we define a temporal energy functional £i{t) and the corresponding dissi¬ 
pation rate 'Di (t) by 


£i{t):=\\ftr+\\ptr + \\u,r+ ^ 

l<|a|<3 ^ 




l + p 


V 


+ ||a“ntf'|>+Kifo(i)+«:2 V [ • a“Vpt dx, 

J |a|< 2 -^«^ 

(t) := \\y{aubuPuUt)rH2 + ||6t -«*||^3 + ^ {||a“{I - P}/t||2 + ||5“Vwif}, 


l«l<3 


where ki, K 2 > 0 are small constants to be chosen later. Obviously, we have 

£ 1 ( 1 ) ~ Il/i(^)lli 2 (ff|) + Il(pt,'*^t)(0llff3- 

According to p.34|) . p.39|) . (I3.40|) . and by choosing ki, K 2 ande snfficiently small, 
we finally obtain that 

piit) + \ E ii5“{i-p}/tii^ 

l«l<3 

+ ^(ll^(Ot,^t,Pt,'*^t)llll 2 + \\ut — btWjfS + ||VUt||^3^ 

<C'|||p||/i-3 -I- ||p||^3 -f ||M||/i-4 -I- ||m||^ 3 -b ||p||/i-3||M||/i-4 -b ||p||^3||u||^3| 

x{ E ||9“{I-P}/t||^ + ||V(at,6t,pt,iXt)||2,2 + ||irt-6*11^2} 

|a|<3 

+ C(l + ||p||^4 + ||ir||^3){ E P“{I-P}/f. + l|V(a,6)11^2 

|a|<3 

+ Il^p|lff3 + II^^IIh 4 + 11^ ~ ^llH3|(||/t||i2(/J3) -b ||(pt, mOIIhs)- 
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For simplicity, by the definition of £i{t),'Di{t) in (I1.14|1 . (I1.15L we can rewrite the 
above inequality as 

^fi(t) + APi(t) < C{£ht) + + C{1 + (3.41) 

Next, we need to estimate the mixed space-velocity derivatives of ft- Since 
||3|P/t|| < C\\d°‘ft\\ for any a,/3, we only estimate ||9|{I - P}/t|l. From (14.231) . 
we easily deduce that 

dt{l - P}ft + V ■ V,{I - P}ft + u ■ V4I - P}ft - ■ u{I - P}ft 

= C{I - P}ft + P{v V,{I - P}ft + u ■ V„{I - P}ft - ■ u{I - P}ft) 

-{I-P}(u-V,P/i+u-V„P/t- iu-uP/t) +{I-P}Gt, (3.42) 

where {I — P}Gt is defined by 
{I-P}G* 

= pC{I - P}ft + ptC{l - P}f - - P}ft + ^pu ■ u{I - P}ft 

- {ptu + (1 + p)ut) ■ (V„ - ^) {I - P}/ 

+ P{pu • (V„ - ^) {I - P}ft + {ptu + (1 + p)ut) ■ (V. - ^) {I - P}/} 

- {I - P}{pu ■ (V. - ^)pft + {ptu + (1 + p)ut) ■ (V„ - ^)P/}. 

Similarly to the proof of Proposition 14.51 we can prove that 

Proposition 3.5. Assume that {f,p,u) is the solution obtained in Theorem \l.l[ 
Then we have 

A ^ iia|{i-p}/tf+ A ^ Wp{i-p}ft\\l 

|/3|=fe \P\=k 

|a| + |/3|<3 |a|-t|/3|<3 

<g{ ^ W'{i-P}ft\\l + \\vht\\]j,-,] + Cx 2 <k<A Y. ll9|'{i-P}/tll" 

|a'|<4-fe l<l/3'|<fe-l 

|a'| + |/3'|<3 

+ c(i + ||p||1^3)||Vu||^J ^ ||a“'{i-p}/tf + ^ ||a^;{i-p}/,f| 

l«'|-F/3'l<3 

+ c(i + Ml. + ||u||?,3){ ^ ||a|;{i - P}f\\l + \\Vb\\l,]\\{p,, 

|a'| + |/3'|<3 

+ Cmm + Mj,. + \\p\\U\nrm){ Y lia“'{I-P}/i||^ 

|a'|<3-fc 

+ ^ \\dp{i-p}ft\\i+\\vhrH.}. 

1<|/3'|<3 
|a'| + |/3'|<3 


(3.43) 
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Multiplying (j3.43|) by suitable constants Rk and taking summation over fc, we 
obtain the following inequality: 

^ ^ ^ \\d^{i-F}ftr+x 

l<k<3 \p\=k 1<I/3|<3 

|a| + |/3|<3 |a| + l/3|<3 

<C{J2 P“{I-P}/tll" + l|V6*f^.} 

|a|<3 

+c{i + \\p\\u\\^n\\i.{j2\\d‘^{i-p}ftr+ iiw-p}/tf} 

l«l<2 1<I/3|<3 

|a| + |/3|<3 

l«'| + |/3'|<3 

+ C(||pf^3+||u||?,3 + ||p||l,3||u||?,3){ ^ ||a“{I-P}/*f, 

|a|<3 

+ ^ \\d^{l-P}M\l + \\VhrH.]. (3.44) 

1<I/3|<3 
|a| + |/3|<3 

Now, we define £2 (t) and £>2 (t) as 

£2{t):= E E IIW-P}/tf, 

l<fe<3 \p\=k 

|a| + |/3|<3 

V2it):= E IIW-P}/tf.- 

1<I/3|<3 
|a| + |/3|<3 

Therefore, according to the definition of 'D 2 {t) in (I1.17L (13.441) can be rewrited as 

+ Ap2(t) < C'Di{t)+C{£i{t)+£f{t)){'Di{t)+t>2{t)) 

+ C(l + £i{t)^ + ^2{£))- (3.45) 

Thus, we define a total energy functional £{t) and the corresponding dissipation 
rate V{t) by 

£{t) := £i{t) + K3f2(t), 

2?(t) := 'Di{t) + K3f>2{,t), 

where K 3 > 0 is a small constant to be chosen later. 

With the aid of (I1.18I) . (|1.19L(I3.41I) and (I3.45p . we conclude that 

^£{t) + \V< C{£^ (t) + £^{t))V + C(1 + £2'D£- (3.46) 

By the uniform a priori estimates obtained in Section 4, we deduce that £(t) is 
sufficiently small. Therefore, we obtain that 

^£{t) + XD< C{1 + £‘^)V£. 


(3.47) 
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Applying Grownwall’s inequality to (I3.47[) . we obtain 

p+oa 

sup £{t) -\- \ / Vdt < CE{Q) < + 00 . 

0<i<C5O JO 

Meanwhile, according to equations (ll.7p - ()1.10|l and the uniform a priori estimates 
obtained above, we deduce that, 

^+oo 

/ (ll/tf + ||ptf + ||u*f)df<+oo. 


Thus, 


r‘+oo 


/ (Il/t|lw3 +||p*||^3 + ||u*||^3)dt<+00. 

Jo 

On the other hand, we can apply Lemma l3.II to obtain 

+ Il“‘llff3 < YTt’ 

Now, we define the following functionals 

|a|<2 ij 

- E / ■ hdx, 

|a|<2"' 


(3.48) 


\/J(l + P) 


1 + p 


d°‘p 


||9“ 


£i{t):=\\fr+\\pr+\\ur+ e ii^“/ii' + 

l<|a|<3 ^ 

+ Ti£o{t)+T2 E /9“u-9“Vpdx, 

|a|<2"' 

Pi(t) := \\W{a,b,p,u)rH. + \\b-u\\j,s+ E {l|5“{I-P}/ll" + l|a“Vuf}, 

|a|<3 

f2(t):= E E IIW-P}/f, 

l<fc<3 |^|:=fe 

|a| + |/3|<3 

f>2{t):= E I|5|{I-P}/|1^, 

1<|/3|<3 
l«l + l/3|<3 

£{t) := £i{t) + Tz£2{t), 

'D{t) ■= f>i{t) + T3V2{t). 

According to the results obtained in Section 4, we finally obtain that 


{ E + ll^ll 


sup 

■ |a| + |/3|<l 

< E Iiv5|/f+ Iivp||?i2 + ||vu||2,. <p(t)<-AY(t) 

|c| + |/3|<2 

<C{\\f\\Hi^ + ll(p,M)||_f/3)(||/t||//3^^ + \\{pt,ut)\\H!>) < c{l +t)“5. (3.49) 
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This completes the proof of Theorem 11.11 

3.2. The case of periodic domain. In this subsection we study the large time 
behavior of classical solutions when 17 is a spatial periodic domain T^. We shall use 
the uniform a priori estimates obtained in Section 4 below. It follows from (14.331) 
together with (14.351) that 

^£{t) + Writ) < C{£^ {t) + £^{t))Vj{t). 
at 

Using the fact that £{t) is small enough and uniformly in time, and £{t) < C'Dj{t), 
we then obtain 

^£(t) + X£(t) < 0 
at 

for all t > 0. This gives the desired exponential decay by applying Gronwall’s 
inequality, and hence completes the proof of Theorem 11.21 


4. A PRIORI ESTIMATES OF THE CLASSICAL SOLUTIONS 

In this section, we shall establish the uniform-in-time a priori estimates in the 
spatial domain 17 = or which have been used in Sections 2 and 3. We need 
the following two assumptions: 

(1) (/, p, u) is the smooth solution to the Cauchy problem (jl.7|) - (ll.l0|) for 0 < 7 < T 
with a fixed T > 0; 

( 2 ) (/, p, u) satisfies 

sup {\\f{t)\\Hi, + Up,u){t)\\H^} <S, (4.1) 

0<t<T 

where 0 < <5 < 1 is a sufficiently small generic constant. 

First, we introduce a lemma which is useful in the subsequent estimates: 

Lemma 4.1 (see m)- There exists a positive constant C, such that for any f,g€ 
and any multi-index 7 with 1 < I 7 I < 4, 

||/||l»( 0 ) < (4.2) 

< C'll/llH 2 (n)l|Vxg||H 2 (n), (4.3) 

\\92if9)\\L^{n) < C\\X/xf\\H3(n)\\Vo,g\\H3{Q)- (4.4) 

4.1. The case of the whole space. In this subsection we deal with the uniform- 
in-time a priori estimates in the whole space 17 = 

Proposition 4.1. For smooth solutions of the problem (I1.7I) - (I1.10I) . we have 

P. + A(||{I - P}f\\l + I|6 - uf + llVirf) 

< C(||(p, u)||^= + \\p\\h3 ||u||^.)(|| V,(«, b, p, u)r + ||u - + ||{I - P}f\\l) (4.5) 


for all 0 < t < T. 
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Proof. Multiplying (ITTll-lfm by /, p, and u respectively and then taking integra¬ 
tion and summation, we get 

+ llpf + l|wf) + j (--^{1 - P}/: /) da; -b J da; -b ||6 - uf 

= J u(^—vf,f'^dx — J a\u\^ dx — J pdivudx — J p'{l)yp-udx 
— J {u-S/u) ■ udx — J V ^ ^ Vm ■ udx — j — p'(l)^ Vp • uda; 

—-y p^divuda;-b JJ p(^Cf — u ■ V yf + -u ■ vf + u ■ v'/M^ f dxdv. (4.6) 
By (I1.12|) . we have 

Thus, we only need to estimate the terms on the right hand side of the equality 
(TO . For the first two terms, by taking the same computation as that in m , we 
get 

J u(^^vf,f'jdx- J a\u\'^dx 

< C(|| + Mh^I - Pf}\\l + C||ix||^.(|| V,(a, &)f + Ilia - fcf). 

Without loss of generality, we can assume that p'{l) = 1, then, one has 

J pdivudx-b J p'{l)Vp ■ udx = 0. 

For the next four terms, using Holder’s, Sobolev’s inequalities and Lemma we 
have 


J{u ■ Vit) • udx < Cj 

V—-—Vu • udx < C| 
1 -b p 

+ -P'(l))vp-udx < C\ 

<C\ 

^ J divudx < C\ 
<C\ 


u|U3||V,u||i2||u|U6 <C||u||^i||Vu||i., 
^llaa^dlVpf-bllVuf), 

P\\L3\\yxP\\L^\\u\\L<^ 

Pllff^dlVpf+ ||Vuf), 

p|lL3||Va,u||L2||p||i6 

PlIiL^dlVpf + ||Vuf). 


Here, we have used the facts that ||u||i6 < C'||Vu||i 2 and ||u||i3 < CHuHj^i. 
Using the macro-micro decomposition (ll.llL we rewrite the last term as 


— u ■ Vyf -b —u ■ vf + u ■ vy M ) / dxdu 


//K« 

= J p(—£{I — P}/, /) dx-b J p{u — b)-bdx+- JJ pu-vp dxdu. 
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It is easy to see that 

J p{u-b)-bdx< C\\p\\L3\\u-b\\L2\\b\\L6 < C'||p||ffi(||V 6|||2 + ||m- 6 III 2 ). 

Noticing that 

(\vj, f) = ab, + {v.Pf, {I - P}/) + (in., |{I - P}/p), (* = 1 , 2 , 3), 

the term ^ JJ pu ■ vp dxdz; can be estimated as follows: 

J pu-badx < C'||p|li 3 ||M||i 3 ||a||i 6 || 6 ||i 6 

<C\\p\\H4u\\mi\\yar + \m^), 


I p{v.Pf, {I - P}/) dx < C\\puU3\\{a, &)l|La|l{I - P}/|| 

<C\\p\\M\H3i\\ViaP)\\^ + \\{l-P}f\\l), 

J |{I - P}/p) dx < C\\pu\\l^ ||{I - P}f\\l 

<C\\p\\H4u\\H4{l-P}f\\l. 


Therefore, the last term is bounded by 


I p{-C{l-P}fJ)dx 

+ C{\\p\\h3 + \\p\\h4u\\h3) (II V(a, 6 )f + ||ir - + ||{I - P}f\\l). 


Plugging all the above estimates into (14.611 and using (I4.1I1 . we obtain (14.51) . 


□ 


Proposition 4.2. For smooth solutions of the problem (I1.7I) - (I1.10I) . we have 


l_d 

2 dt 


l<|a|<4 


1 + P 



+ A E {l|{I-P}5“/ll^ + l|9“(6-u)f+ ||V9“«|p} 

l<|a|<4 

< C'{h||^4 + ||p||^4||zi||^4(l + IIHI^,) + ||H|i44(l + IIHII,.)} 

x{\\VPa,b,p,u)\\%s+ E ll{I-P}5“'/ll"} 

1<I“'I<4 


(4.7) 


for all 0 < t < T. 

Proof. Applying differentiation 9“(1 < |a| < 4) to the system (ll.7l) - (ll.9|) . we have 
at(5“/) + V ■ v,(a“/) + u ■ v,(a“/) - d^u ■ vVm - zd^f 

= + p)u ■ vf\ + [-9“, u ■ V,]/ + 9“{p(£/ -u-wp + u- vPm)}, 

(4.8) 

9t(9“p) + u- Vi9“p+ (1 +p)divcl“it = [-8^^, pV x-]u + [-d°‘,u- Vx\p, (4.9) 

5i(9“u) + u • V((9“u) + - u) 

1 + p Vl + P / 
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= V,]m + 


-d° 


p'(i + p) 

1 + p 



p — d‘^{ua), 


(4.10) 


where [A, B] denotes the commutator AB — BA for two operators A and B. Now, 
multiplying (I4.8|) - (l4.10p by 9“/, and d°‘u respectively and then taking 

integration and summation, we have 




\/p'(l + P) 


1+ P 


d°‘p 




ic^.,l|2 


+ a“/) dx + y J^\\7{d°^u)\'^dx+\\d°‘{b 

= y([-5“,u-v.]/,a“/)dx + yi(a“[(i + p)u.u/],a“/)dx 


+ J[-d‘^,P'^a:-]uj^^j^^d°‘pdx + Jd’^pdx 
+ A—9“, u • Va,]Mi9“Mda; — i f |cl“itp divuda; 


+ 



P'i^ + P) 

1 + p 


Vx]pd°‘udx 


8°" pd°‘u 


^ p'[l + p) 

1 + p 


dx 


1 

2 


|9“p|2div ( 


y(l + p) 

(1 + P)^ 



d°‘{ua)d°‘udx 


1=1-^ l<P<a 

+ y (a“{p(/:/ -u-y,f + u- vVm)}, d^f) ^^ + \j d^ 

14 

:=^/„ (4.11) 

1=1 


where Ca,p are constants depending only on a and /3. Each term in (14.111) can be 
estimated as follows. For /i,/ 2 ,/ 5 , and Iio, we can carry out similar arguments to 
the proof of Lemma 2.3 in [T^] to obtain (the details is omitted here) 


h <c||Vu||ff3||v,/iU2(ff3)||v.a“/||, 

/2<C(1 + ||p||^4)||Vu||^,3||V,/|U3(^3)M“/||, 
h<C\\Vu\\]j,Wul 
Iw<C\\Vu\\HA\^a\\mWu\\. 

Using Holder’s, Sobolev’s, and Young’s inequalities, we easily get the following 
bounds: 

h + h< c||Vu||j^3||Vp||^3||a“p|| < c||vu||j^3||vp||^3, 
h < qi divM||ioe||a“Mf < c||vu||ff3||a“uf, 

/7<C'(||Vp||^3+||Vp|||,3)||Vw||^3, 

h<c\\u\\HA\yp\\l.. 
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I 9 < Cdl divM||L=o + ||M||L~||Vp||L=)||a“pf 
<c(i + ||p||^3)|H|^3||a>f, 

hi < C||Vp||ff2||va“w||||9%|| < Cd|Vp||^3||5“«||2 + e||V9“uf, 

/12 < C||Vp||ff3||V9“w||||9%|| + C(||VHI//3 + ||Vpr^3)||Vwf^3 

< C,(||Vp||ff3 + ||VH||^3 )||Vw|| 2,3 +e||V9“^zf 

with e > 0 a small constant. 

By means of (113 and (03, one has 

sup \dtp{t,x)\ < (1 + ||/o||l=°)|| divu||L=o + ||m||l=o||Vp||l=o < C\\u\\h3. 

0<t<T,a:eR3 

Then, the term /14 can be controlled by the following bound 


Now we estimate the term / 13 . We have 
JJ-d°‘{pu-yyf)d°‘fdxdv = jj d°‘{puf)-V^d^fdxdv 

<c||vp||,43||vu||,^3||v,/iu2(,,3)||v,a“/||, 

Jj 5“(H • v\/Md°‘fdxdv < C\\d°‘{pu)\\\\v\fMd°‘f\\ 

<C\\^p\\H3\\^u\\H3Wf\l. 


JJ d°‘{pCf)d°‘fdxdv 


d°‘{pg)d°‘gdxdv 


<C||Vp||^,3||Vg||2,(^3) 
<C||Vp||,^3 ^ ||5“7||^ 


l<|a'|<4 

with g = •\/MV^(M“U)- Thus, /i 3 can be bounded by 

C{\\p\\H^ + \\p\\H4u\\H^){\\V,{aAu)rm+ ll{I-P}5“'/ll^}. 

l<|a'|<4 

Plugging the estimates on /i (1 < i < 14) into (14.1111 and taking summation over 
1 < < 4, we obtain 03. □ 


In order to get the energy dissipation rate || V 2 ;(a, 6)||j:/3, we need to study the 
equations satisfied by a and b. We follow some ideas developed in [T2] where the 
lower order estimates similar to Proposition 14.31 below were obtained. It is easy to 
verify that a and b satisfy the following equations: 

5ta + div6 = 0, (4.12) 

dtbi + dxid + — P}/) = —(1 + p)bi + (1 + p)ui{l + a), (4-13) 

j 

dxjbi + dxfbj - (1 + p){uibj + Ujbi) = - P}/ + Tij{l + r + s) (4.14) 
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for 1 < i,j < 3, where Fij is the moment functional defined by 

^i,j = {{viVj - 1)'/M,g) 
for any g = g{v), and s are defined respectively by 
/:=-i;-V4I-P}/ + £{I-P}/, 

r-.= -u- V,{I - P}/ + • «{! - P}/, 

s := pM-^Vy ■ - P}/ + VM{1 - P}/ - uy/M{l - P}/). 

In fact, (I4.12|) and (14.131) can be obtained straightforwardly by multiplying (11.71) by 
'/M and Vi^/M (1 < i < 3) respectively and then taking the velocity integration 
over . To obtain (I4.14|) , we can rewrite (11.71) as 

dtPf + V ■ V,P/ + u ■ V„P/ - in • vPf + Pi/ = -dt{l -P}f + l + r + s, 

then apply Fy to it and use (j4.12l) . 

We have the following estimate. 

Proposition 4.3. For smooth solutions of the problem (I1.7I) - (I1.10I) . we have 
Afo(t) + A||V,(a, b)\\jj, < C{\\{I - P}/lli.(ff4) + ||n - b\\l,) 

+ C[\\p, u||^3 + llpllffsllnll^s) 

X (||V,{I - P}/Ili.(i^3) + ||n - b\\],s + II V,(a, b)fH,) 

(4.15) 

for all 0 < t < T. 

Proof. It follows from ()4.14l) that 

X a“[(l + p){uibj + Ujbi) - 9tFij({I - P}/) + rij{l +r + s)] da; 

= - ^ E / - P}f) dx 

+ Y. [ + djdth)d‘^r,,j{{l - P}f) dx 

+ ^ y* d'^idibj + djbi)d'^[{l + p){uibj + Ujb^) + Tij{l + r + s)] dx. (4.16) 
Applying (I4.13|) , Lemma 14.11 and Young’s inequality, we obtain 

X] J d^^idAb, + - P}/) dx 




<e||V.a||l,3+C.||V4I-P}/||i.(^3) 
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+ C'(l + ||p|||f3)(||M — bW^S + ||'«|||f2||Va;a|||f2) 

with e > 0 sufficiently small. For the final term on the right hand side of (I4.16L 
we have the following estimate: 

f d‘^{dibj + djbt)d°'[{l + p){uibj + Ujbi) + + r + s)] da; 

2J' 

+ C' + P){'u-ibj + Ujbi)\\‘^ + \\d°'Tij{l + r + s)||^). 

i ,3 

According to Lemma |4.1l the definition of Tij, and the expressions of I and r, we 
get 


E < C{1 + ||p||?^a)||u||?^3||V,6||2,3, 

Ell^“Ua(0f <OI{I-P}/lli3(^4), 

E l|5“Ua(Of < C||«||2,3||V4I - P}/Ili2(^3). 

For r2^j(5), we have the following estimate: 

rij('S) ^ ~P J + (Vidvj + v,dy^) - {uiVj + UjVi)}'/M{l - P}/dv. (4.17) 
Now, we deal with the terms on the right hand side of (I4.17p . First, we have 


a“ J pUjV^y/M{l - P}fdi 


da; 


J - P}/)da 

<C j J \d‘^{pu,{l-P}f)\^dxdv 
<0|P|I?^3||OI?^3||V4I-P}/Ili2(^3). 

Similarly we obtain that 

j pu,u,VM{I-P}/d^;||Eq|p||?^3||w||l,3||V4I-P}/|U2(ff3), 

5“ I pv.v,VM{ 1 - P}fdvf < Cllpll^all V.{I - P}/lli2(^,3), 

J +u.a,,)fM{i-p}/du||Eq|p||?f3||v4i-p}/||23(^3). 


E \\9‘^i^^bJ + d,bi)\\^ = 2|| V,5“6f + 2|| V, • 9“6f. 

-ij 


Notice that 
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Using the above equality and plugging the above estimates into (I4.16|l . and then 
taking summation over |q;| < 3, one gets 

/5“(a,&, + 9,6,)a“r.,,({i-P}/)dx + ||v,5“6f+ ||v,-a“6f 

™ |c<|<3 ^,j 

< e||Va||^3 + Cellll ” P}/|||2(/j4) + C\\u — &||^3 + C{\\{p, u)||^3 + IIpII^s|| u||^3) 
X (||V,{I-P}/||i2(^3) + ||u-&||l,3 + ||V,(a, 6)11^3). (4.18) 

On the other hand, by means of (I4.12p and (j4.13|) . it follows that 

||a“Vaf = -A f 5“a5“div6da;+|ia“div6f 
dt J 

+ ^ f d°‘diad°‘^{l + p){ui - h) - - P}/ + (1 + p)uia| dx 

* 3 

^“«^“div&dx+||a“div6f+ i||a“Vaf+ C||V,{I-P}/||i2(^3) 


■ C'(l + I|p|lw3){ll'“ ~ + ll“llfI3||Va||^3}. 


(4.19) 


_d 

dt 


Summing (14.191) over |a| < 3 and taking e = -j, and then adding the results into 
(14.181) implies (14.151) . □ 

Proposition 4.4. For smooth solutions of the problem (I1.7l) - (ll.l0p . we have 

• 9“Vpdx + A||Vp||^3 

|c.|<3-^«= 

< C(||u - 611^3 + ||Vu||^4) + C(||p||^4 + 11^111,4 + hl|?,3)||V(p,u)||^3 (4.20) 

for all 0 < t < T. 

Proof. Taking differentiation 9“ (|q;| < 3) to (II.9L and by carrying an direct calcu¬ 
lation, we get 

p'(l)||va“pf =- J yd°‘pd°‘dtudx + J Vd°‘pd°‘{b-u)dx 


+ J - u • Vu -t- - 

■= Qi + Q 2 + Qs- 


p'J3+_p) 

i + p 


-p 


'(i)]vp} 


dec 
(4.21) 


For Qi {i = 1, 2,3), applying (II.8L Holder’s, Sobolev’s and Young’s inequalities, we 
have 

Qi = - ^ y Vd°‘pd°‘u dx -f y divu5“[(l + p)dWu + u- Vp] dx 
< -^y Vcl“pi9“udx-I-0119“ divulp-I-C|| p||//4||Vu||^3, 
Q2<^||va“pf+ c||u-&f^3, 

Qs < ^11 Va“pf + C||u • Vu||^3 + C 


1 + p 


-Au 


ff3 
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< 


+ c 

y(l) 


p'il + p) 


l + p 




Vp 


wvayr + c||u||^3||Vur^3 + c\\p\\H4vdyr + c^livu|| 


-C{\\p\\ 




I?I3)||(VP,VU)||^3. 


2 

H* 


With the help of (14.IL plugging the above estimates into (I4.21L we obtain (14.201) . 

□ 


Reorganizing the estimates obtained above in the Propositions 14.Hl4.4l one gets 

^£i(t) + XVi(t) 

< C{\\{p,u)\\h^ + ||p||/i-4||u||/i-4(l + |1 p||^ 4 + ||p||/i-4||u||/i-4) + ||p||^4(l + ||p||^4)} 

X |||V(a,6,p,u)||^3 + ||6-u||?,4+ ^ ||a“{I-P}/||2|. (4.22) 

^ l«l<4 ^ 

Next, we need to estimate the mixed space-velocity derivatives of /, i.e., d^f. 
Since ||9|P/|| < CIlS"/!! for any a and /3, we only need to estimate ||9|{I — P}/|| 
below. Let us apply I — P to both sides of (13 to get 

dt{l -P}f + v- V,{I -P}f + u- V„{I - P}/ - iu • u{I - P}/ 

= £{1 - P}/ -f P{u • V,{I -P}f + u- V„{I - P}/ - iu • u{I - P}/} 

- {I - P}{u • V,P/ + u ■ V„P/ - iu • vPf} + {I - P}G. (4.23) 

where {I — P}G is defined by 

{I - P}G := p[c{l - P}f + ■ u{I -P}f-u- V4I - P}/ 

+ P(u • V„{I - P}/ - iu • u{I - P}/) 

-{I-P}(u-V,P/-iu-uP/)}. 

Now we give a simple derivation of the equality (14.231) . First, 

{I - P}(u • V,/) = V ■ V,{I -P}f + vPf- P{v ■ V,{I - P}/) - P(u • V,P/) 
= V ■ V,{I - P}/ + {I - P}(u • V,P/) - P(u • V,{I - P}/). 

We can deal with the terms {I — P}(u • V„/), {I — P}(u • vf) in the same way. 
Meanwhile, by the definitions of £andP, one has 

{I-P}{u■v^/M)=0, {I-P}Cf = C{l-P}f. 

Then (14.2311 follows. 

In the proof of the following Proposition, we adopt some ideas from [161 Lemma 
4.3]. 
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Proposition 4.5. Let 1 <k < A. For smooth solutions of the problem (ll.7|) - (jl.lOD . 
we have 

A ^ ||a^{i-p}/f+ A ^ \\d^{i-p}f\\l 

\0\=k |/3|=fc 

|a| + |/3|<4 |a| + |/3|<4 

<c(i + \\{p,u)\\j,, + MUHU Il5“'{i- P}/ll^ 

|a'|<4—fc+1 

+ c(||p||^3 + iiiiii^a + MUWWU E \\dp'di-P}f\\l 

1<I/3'I<4 
l«'| + |^'|<4 

+ + C{1 + ||p||l,3)hf^4-.||V6||^3 

+ Cx2<fc<4(l + ||pfff3) E \\di{l-P}f\\l (4.24) 

l<l/3'l<fc-l 
|a'| + l/3'l<4 

for all 0 < t < T. Here xe denotes the characteristic function of the set E. 


Proof. Fix fc (1 < fc < 4). Choosing a and (3 such that |/3| = k and |a| + |/3| < 4, 
multiplying (14.231) by — P}/ and then taking integration, one has 


l_d 

2dt 


||9|{I - P}/f + / {-Ld^{l - P}/, {I - P}/) dx := E 


(4.25) 


2=1 


with 

Ji=j , u • V,]{I - P}/, d^{l - P}/) dx, 
j2 = J {d:[d^,-\vm - p}/, d^{i - p}/) dx, 

J 3 = j {-d${u ■ V„{I - P}/), a|{I - P}/) dx, 

Ji = J • HI - P}/), 5^1 - P}/) dx, 

J5 = j (5|p(u-V,{I-P}/ + u-V4I-P}/ 

-iu-u{i-p}/),a|{i-p}/)dx, 

J6 = J {-d^{I-P}(v V,P/ + u • V,P/ - iu • uP/) ,d^{I- P}/) dx, 
Jr = J{d^{l-P}G,d^{l-P}f)dx. 


Here the fact that [d^,C] = [d^,—\v\'^] has been used. 

Now we start estimating each term in (14.251) . For the terms (i = 1, • • • ,6), 
we have 


Ji < p\\d^{i- p}/f + V • v„]a“{i- p}/f 
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<, 7 l|a|{i-p}/f+ c, ^ ||5“'v4i-p}/f 

\a.' I <4—fc 

+ X2<k<ACrj E ii5|'{i-p}/f, 

l<l/3'l<fc-l 
l«'| + l/3'l<4 

J2 < vmi - p}/f+ cj[d^,-\v\^m{i - p}/f 

<i 7 ||a^{I-P}/f+ c, ^ ||5“'{I-P}/||2 

|a^ I <4—fc 

+ X2<k<ACrf 

l«'| + l/3'l<4 

J3 < ^Iia^{i - P}/f + c,||5“(u • v.af {I - p}/)f 
<r;||a^{I-P}/f+ C,||w||^3 E P^''{I-P}/f, 

1<I/3'I<4 
l“'l + l/3'l<4 

J4 < vWpil - P}/f + c,|ia“(« • dl{v{i - p}/))f 

<,7ii9|{i-p}/f+c',hi|2,3 E 

1<I/3'I<4 
l“'l + l/3'l<4 

+ E l|5“'{I-P}/f, 

|q(^|< 4—fc 

Js < r^\mi - P}/f + C,|15|P(i; • V,{I - P}/)f 

+ C,||a|P(u • V4I - P}/)f + Cr^WpViu ■ HI - P}/)f 
<r;||5|{I-P}/f+ C, ||V,a“'{I-P}/f 

|a'|<4-fc 

+ ^ lia“'{i-p}/f, 

|q^ I <4—fc 

H - P}/f + c,||a|{i - PIH • V,p/)f 

+ c,||a|{i - p}H • v.p/)f + c.iiaHi - p}H • «P/)f 

< P}/f + C,(||V6||^4-. + ||V6||^3h||^4-.)- 

For the term J 7 , we have the following calculation and estimates: 

J7 = y'(a|(pHi-P}/),5Hi-P}/)dx 

+ lj{d^ {pu ■v{l-P}f),d^{l- P}/) dx 
- |(a^(pu-vHi-p}/)H^{i-P}/)dx 
+ J (9|p(pyV„{i-P}/-ip4i-Hi-P}/),5|{i-P}/)dx 
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- J (^d^{l-P}(^pu-y.Pf-^pu-vPf),d^{l-P}f)dx 

5 

i=l 

We can adopt the above similar estimates to deal with Yi{2 < i < 5). Thus, we 
only give the following bounds: 

Y2<p\\d^{l-P}ff 

+ cjp\\uwfHJ P“'{i-p}/f + E 

^ |a'|<4-fe 1<|/3'|<4 ^ 

|a'| + |/3'|<4 

Y3<p\\d^{i-p}fr+cjp\\uH%, E iiai'{i-p}/f: 

1<|/3'|<4 
|a'| + |/3'|<4 

Y4<v\\d^{I-P}ff+Cr,\\p\\UH%, E l|5“'{I-P}/f, 

|a'|<4-fe 

>5 < v\\d^{i - P}/f + cJpWI, Ml.., M.bM- 

For Yi, we give a detailed calculation: 

Yi = II d:{pLd^{l - P}f)d^{l - Pjfdxdv 

+ II d^{p[d^,-M]{l-P}f)d^{I-P}fdxdv, 

= II pLd^{i-P}fd^{i-P}fdxdv 
+ E [[ d^pLd^~^{l-P}fd^{l-P}fdxdv 

l< 7 <a 

+ II a“(p[af,-|up]{i-p}/)a^{i-p}/dcrdu, 

:=rii+ri 2 + ri3. 

For Yii, we can move it to the left hand side of the equality (I4.25F Thus, we only 
need to deal with Y 12 and Y 13 . We have 

>"12 = - II d'^pVMV,{M-^d;-^{I-P}f)VMV,{d^{l-P}f) dxdu 
<C||a|{I-P}/||.||p||^3 E \\dp'{i-'P}fh 

|/3'|=fe 
l«'l + l/3'l<4 

<C\\p\\h^ E \\dp'>{'^-P}f\\l 

\p'\=k 
|a'| + |/3'|<4 

^13 < - P}/f + cj[d^, -kP]a“p{i - p}/f 
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<,7l|5^{i-p}/f+ c,ihi^ 3{ ^ ||a“'{i-p}/ii2 

\ol' I <4—A: 

+X2<fe<4 ^ \\d^:{i-p}f\\i}. 

|a'| + |/3'|<4 

For the second term in the left hand side of the equality (14.251) . one gets 

J{1 + p){-Ld^{l - P}/, 9|{I - P}/) dx 

> Ai||{I-Po}/||^ 

> - P}f\\l - Ai||Poa|{I - P}f\\l 
>^\\d^{l-P}f\\l-C\\d^{l-P}ff. 

Plugging all the above estimates into (14.2511 and choosing r] sufficiently small, we 
obtain (14.2411 . □ 

Remark 4.1. According to the above Lemma, we can choose some suitable con¬ 
stants Cfc, such that 

4 E E ii5|{i-p}/f+A E \\dp{i-p}f\\i 

^ l<fe<4 |/3|=fc 1<|/3|<4 

l“l + l^l<4 |a|-l-|^l<4 

<c{\\{p,u)\\j,, + \\p\\UHU E l|a“{i-P}/f. 

|a|<4 

+ C{\\p\\h^ + \\u\\1. + \\p\\UW\\U E \\d^{i-P}f\\l 

1<II/3|<4 
a\ + \p\<4 

+ C'hlll4a(l + ||p||l,3)||V6||?,a+c(||V6||^3+ E I|5“{I-P}/I|4- (4-26) 

l«l<4 


With the aid of the inequalities (I4.22L (14.261) and (14.11) . we have 

+ XV{t) < C{£i(t) + £^{t))V{t) < C(S + S^)V{t). (4.27) 

So, as long as 0 < i5 < 1 is sufficiently small, the integration in time of ()4.27p yields 

£{t) + X [ V{s)ds < £{Q) (4.28) 

Jo 

for all 0 < t < r. Besides, (I4.1|l can be justified by choosing 

^(0) ^ + ll(Po,uo)||^4 


sufficiently small. 
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4.2. The case of periodic domain. In this subsection we deal with the uniform 
a priori estimate when II is a spatial periodic domain T^. Using the following 
conservation laws in the case of torus, 


_d 

dt 


Fda;du = 0, — / nda: = 0, 

dtj 


— I nudx + 


vF dxdv J “ 

and by the assumptions of Theorem 11.21 it follows that 

J a dx = 0 , J pdx = 0 , j (6 + (1 + p)m) da; = 0 


(4.29) 


for all t > 0. 

Thanks to the Poincare inequality and the conservation laws (I4.29|l . we have 


|la|U2 <C||Va|U2, |lp||i2 <C||Vp|U2, (4.30) 

\\u + b\\L 2 < \\b + u + pu\\l 2 + \\pu\\l 2 

< C\\V [b + u + pu)\\l 2 + \\u\\l’^\\p\\l 2 

< C||V(6,u)|U2 + C\\u\\h4^p\\l^ + C\\p\\h4^u\\l^. (4.31) 


Similarly to the whole space case, we can obtain the following estimates: 

^£:i(t) + XVi{t) < C{£^^ +£!)'D4t), (4.32) 

^£2{t) + XV2(t) < CViit) + C{£i + £l)Vi{t) + C{£^ + £l)V2{t). (4.33) 

According to the definition of X>t,i(0i have 

25T.i(i)~ E ll{I-P}a“/||^ + |l(a,&,p,u)||?,4. (4.34) 

|a|<4 

Combining (I4.30|) . (14.311) and (14.321) together, we conclude that 

^£i{t) + APT,i(t) < C{£^ + £!)Vr,4t). (4.35) 
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